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We extend Schwinger's ideas regarding sonoluminescence 
by explicitly calculating the Bogolubov coefficients relating 
the QED vacuum states associated with changes in a dielec- 
tric bubble. Sudden (non-adiabatic) changes in the refractive 
index lead to an efficient production of real photons with a 
broadband spectrum, and a high-frequency cutoff that arises 
from the asymptotic behaviour of the dielectric constant. 

PACS: 12.20.Ds;77.22.Ch; 78.60. Mq 



Introduction: Sonoluminescence occurs when acoustic 
energy induces the collapse of small bubbles, and the 
collapse of these bubbles results in a brief intense flash 
of visible light Qj. There are several competing mecha- 
nisms proposed to explain this phenomenal concentration 
of kiloHertz acoustic energy into PetaHertz electromag- 
netic energy. An interesting mechanism, originally pro- 
posed by Schwinger is based on changes in the zero 
point fluctuations of the QED vacuum. In his model, 
Schwinger estimates the static Casimir energy of an ex- 
panded dielectric bubble, compares it with the Casimir 
energy of a collapsed dielectric bubble, and argues that 
this Casimir energy difference (difference in zero point 
energies) would be converted into real photons during 
collapse of the bubble. This model is often described in 
terms of the dynamic Casimir effect (even if the calcu- 
lation is quasi-static). Several authors have argued that 
this model is not relevant to SL (cf. || and references 
therein). Nevertheless, we feel that the underlying idea 
of SL as a QED vacuum effect has been prematurely dis- 
carded without sufficient analysis. 

In this Letter we consider a variant of Schwinger's 
proposal which is obtained by focusing our attention on 
changes in the refractive index rather than on the bub- 
ble motion. We explicitly compute the Bogolubov coef- 
ficients relating two vacuum states characterized by two 
different values of the refractive index. Non-trivial Bo- 
golubov coefficients imply the production of real photons. 
The spectrum is qualitatively compatible with those ex- 
perimentally observed. Calculations are most easily car- 
ried out for extremely large bubbles (large compared to 
the cutoff wavelength), where the Bogolubov coefficients 



take on particularly simple forms in terms of delta func- 
tions. The spectrum and total energy emission are an- 
alytically calculable. For finite bubbles, the delta func- 
tions are smeared by finite-volume effects and the spec- 
trum can be written down as an integral over a suitable 
sum of spherical Bessel functions. This integral must 
be evaluated numerically and can then be compared to 
both the large-volume estimate, and to the experimen- 
tal situation. Even with a rather crude (step function) 
model for the refractive index as a function of frequency 
the resemblance between observed and predicted spectra 
is quite reasonable. For technical details of the compu- 
tation, and additional discussion of the history of this 
proposal, see ||. 

Basic Features: One of the key aspects of photon pro- 
duction by a space-dependent and time-dependent refrac- 
tive index is that for a change occurring on a timescale r, 
and a photon of frequency u>, then in the high frequency 
limit the amount of photon production is exponentially 
suppressed by an amount exp(— ujt). The adiabatic ap- 
proximation always implies such a suppression j3|. The 
importance for SL is that the experimental spectrum is 
not exponentially suppressed at least out to the far ul- 
traviolet. Therefore any mechanism of Casimir-induccd 
photon production based on an adiabatic approximation 
is destined to failure: Since the exponential suppression 
is not visible out to w « 10 15 Hz, it follows that if SL 
is to be attributed to photon production from a time- 
dependent dielectric bubble {i.e., the dynamical Casimir 
effect) , then the timescale for change in the dielectric bub- 
ble must be of order a femtosecond. Thus any Casimir- 
based model has to take into account that it is no longer 
the collapse from i? max to i? m i n that is important. 

The SL flash is known to occur at or shortly after 
the point of maximum compression. The light flash is 
emitted when the bubble is at or near minimum radius 
-Rmin ~ 0.5 ^m = 500 nm. Note that to get an order 
femtosecond change in refractive index over a distance of 
about 500 nm, the change in refractive index has to prop- 
agate at relativistic speeds. To achieve this, we must ad- 
just basic aspects of the model: We will move away from 
the original Schwinger suggestion and we will postulate 
a rapid (order femtosecond) change in refractive index 
of the gas bubble when it hits the van der Waals hard 
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core. The underlying idea is that there is some physical 
process that gives rise to a sudden change of the refrac- 
tive index inside the bubble when it reaches maximum 
compression. Given the fact that the timescale of such 
a change is much shorter than that typical of the bubble 
collapse we shall consider the bubble radius as fixed and 
equal to the minimal one i? m in- For the sake of simplicity 
we take, as Schwinger did, only the electric part of QED, 
reducing the problem to a scalar electrodynamics. The 
equations of motion are 



\7 2 E = 0. 



(1) 



We shall consider two different asymptotic configurations 
for the gas inside the bubble. An "in" configuration with 
refractive index ?ij n , and an "out" configuration with a 
refractive index n ou t- These two configurations will cor- 
respond to two different bases for the quantization of 
the field. The two bases will be related by Bogolubov 
coefficients in the usual way. Once we determine these 
coefficients we easily get the number of created particles 
per mode and from this the spectrum. Using the inner 
product 



=t / e(r,t) # do <fe d 6 x, (2) 

the Bogolubov coefficients are defined as 

aij = {E?»\Ef), i3 = {E^\Ef). (3) 

We focus on the coefficient because it is related to the 
spectrum, number, and energy: 
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(4) 
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Homogeneous dielectric: In the infinite volume limit 



the eigen-modes are plane waves: 

1 exp(i[k • x — ojt]) 



E(x,t) 



(2tt) 3 / 2 



(7) 



We now introduce a "pseudo-time" parameter by defining 
d/dr = e(t)d/dt, that is, r(t) = J dt/e{t). Then 



dr 2 



E = c 2 e{t)V 2 E. 



(8) 



Now pick a convenient profile for the refractive index 



<t) = |(n 2 n + n 2 out ) + |(n^ ut - nfj tanh(r/r ). (9) 

Here To represents the typical (pseudo-time) timescale for 
the change of the refractive index. After computation 
of the Bogolubov coefficient we have to convert back to 



physical time. Defining (n 2 ) = (nf 



t )/2 we get: 
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(2tt) 3 

sinh 2 (n\nf n uj in - ng ut uw| t /(2(n 2 ))) 
sinh (tt n? n iv in t /(n 2 }) sinh (tt n 2 ont u) out t / (n 2 )) 



(10) 



where £o is now the physical timescale of the change in 
the refractive index. For the particular temporal profile 
we have chosen for analytic tractability this evaluates to 
to = i To (nf n + n 2 ut ^j . The sudden approximation is valid 
provided 
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27rf n out max{n in ,n out }' 



(11) 



Thus the frequency up to which the sudden approxima- 
tion holds is not just the reciprocal of the timescale of the 
change in the refractive index: there is also a strong de- 
pendence on the initial and final values of the refractive 
indices. This implies that we can relax, for some ranges 
of values of n; n and n out , our figure of to ~ O(fs) by up 
to a few orders of magnitude. Unfortunately the precise 
shape of the spectrum is heavily dependent on all the 
experimental parameters (K, ni n , n ou t, R)- This discour- 
ages us from making any sharp statement regarding the 
exact value of the physical timescale required in order to 
fit the data. In the region where the sudden approxima- 
tion holds the various sinh(a;) functions in equation (|Io| ) 
can be replaced by their arguments x. Then 
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In the real physical situation n ln is a function of w; n and 
n ou t is a function of w out . Schwinger's sharp momentum- 
space cutoff for the refractive index is equivalent to the 
choice 



n(fc) = n Q(K - k) + 1 &(k - K). 



(13) 



More complicated models for the cutoff are of course pos- 
sible at the cost of obscuring the analytic properties of 
the model. Taking into account the two photon polariza- 
tions 
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The spectrum is 
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Thus at low frequencies, where the sudden approximation 
holds strictly, the spectrum should show a polynomial be- 
haviour (instead of the linear one expected for a thermal 
distribution) . Q The number of emitted photons is 
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The total emitted energy is approximately 
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(16) 



(17) 



So the average energy per emitted photon is approxi- 
mately (E) = jh w mal ~ 3 eV. To extract some 
numerical estimates, recall that in our new variant 
of Schwinger's model we have #Ught-emitting-region ~ 
R min w 500 nm. We take K ohscrvcd = Kn Vlquid w 
27r/(200 nm) so that iT bserved-R ~ 57r « 15. To get about 
one million photons we now need, for instance, rii n rs 1 
and riout ~ 12, or n ln w 2 x 10 and n out ~ 1, or even 
"out ~ 25 and nj n rs 71, though many other possibili- 
ties could be envisaged. Note that the estimated values 
of n° g H and extremely sensitive to the precise 

choice of cutoff, and the size of the light emitting region. 
More systematically, using K o b S0Ive dR « 15 we get 
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Solving for Tl'm 3S 3 function of n ou t and N, and taking 
N = 10 6 , the result is plotted in figure (|l|). For any 
specified value of n ou t there are exactly two values of n; n 
that lead to one million emitted photons. 




n(out) 



1 This statement is independent of the explicit form of the 
profile for the change of the refractive index. Only very special 
profiles (exponential and never-ending in time) provide an 
exactly thermal spectrum M . 



FIG. 1. The initial refractive index n 1IL plotted as a function 
of n out when one million photons are emitted in the sudden 
approximation. 

Finite volume effects: In finite volume the eigenmodes 
are combinations of Bessel functions and Spherical Har- 
monics, subject to 



El — ebubblo-contcnts if T < R, 
(2 = Emedium li V > R, 



(19) 



and satisfying appropriate junction conditions at r = R. 
We limit ourself to just quoting the key result ||. Intro- 
ducing An = < n as - n™\ 
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Where W(f, g) is the Wronskian and the AJf are calcu- 
lable coefficients depending on the frequency, refractive 
index of the bubble, and that of the ambient medium. 
The above is a general result applicable to any dielec- 
tric sphere that undergoes sudden change in refractive 
index. This expression is far too complex to allow a 
practical analytical resolution of the general case. For 
the specific case of sonoluminescence, we have developed 
suitable numerical approximations. In the infinite vol- 
ume limit there were two continuous branches of values 
for TZg^g and n° a * that led to approximately one million 
emitted photons. If we now place the same values of 
refractive index into the spectrum obtainable from the 
Bogolubov coefficients derived above, numerical integra- 
tion again yields approximately one million photons. The 
total number of photons is changed by at worst a few per- 
cent, while the average photon energy (3/4 times the cut- 
off energy) is almost unaffected. (Some specific sample 
values are reported in Table I.) The basic result is this: 
as expected, finite volume effects do not greatly modify 
the results estimated by using the infinite volume limit. 
Note that hu> max is approximately 4 eV, so that average 
photon energy in this crude model is about 3 eV. 
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(E)/hu max 


2 x 10 4 


1 


1.06 x 10 6 


0.803 


71 


25 


1.00 x 10 6 


0.750 


68 


34 


1.06 x 10 s 


0.751 


9 


25 


0.955 x 10 b 


0.750 


1 


12 


0.98 x 10 B 


0.765 
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Table I: Some typical cases. 



In addition, for the specific case rj£L = 2 x 10 4 



= 1, 



we have calculated and plotted the form of the spectrum. 
We find that the major result of including finite volume 
effects is to smear out the otherwise sharp cutoff coming 
from Schwingcr's step-function model for the refractive 
index. Other choices of refractive index lead to quali- 
tatively similar spectra. These results are in reasonable 
agreement with experimental data. 




FIG. 2. Spectrum dN/dx obtained by integrating the 
approximated Bogolubov coefficient. For n out = 1 and 
R = 500nm the relation between the non-dimensional quan- 
tity x and the frequency v is x ~ v ■ 10.5 -lO" 1 ^. Soisi 11.5 
corresponds to v » 1.1 PHz. The curve with the sharp cutoff 
is the infinite volume approximation. Finite volume effects 
tend to smear out the sharp discontinuity, but do not greatly 
affect the total number of photons emitted. 

Conclusions: We suggest that the key to the SL phe- 
nomenon is not the details of the bubble collapse, but 
rather the way in which the refractive index changes as a 
function of space and time. Sudden changes in the refrac- 
tive index will lead to efficient conversion of zero point 
fluctuations into real photons. Fitting the details of the 
observed SL spectrum then becomes an issue of building 
a robust model for such sudden changes of the refractive 
index of the entrained gases as functions of frequency, 
density, and composition. 

A viable conjecture is that the refractive index of the 
trapped gas undergoes major changes near the moment 
of maximum compression, when the molecules in the gas 
bounce off the van der Waals' hard core. Femtosecond 
timescales for changes in the refractive index have al- 
ready been envisaged in the literature fa], but we feel 
that the present model should encourage a more detailed 
experimental investigation regarding the dynamics of the 
refractive index, both that of the entrained gas as well of 
the water surrounding the bubble. 



We stress that any mechanism that provides sub- 
picosecond timescales for the change of the refractive in- 
dex would imply an important contribution of Casimir 
photons in sonoluminescence. It is not inconceivable that 
for some such mechanisms (e.g., sudden ionization) and 
for some regions of the parameter space the dynamical 
Casimir effect could play an important role in addition to 
or in opposition to other proposed explanations of photon 
emission (e.g., Bremsstrahlung or shock waves). 

Generic features for testing our proposal are j3), non- 
thermal behaviour of the spectrum at low frequencies, 
absence of hard UV photons (no dissociation of the water 
molecules) , emission of photons bounded in angular mo- 
mentum (I < KR). Detection of correlations M in the 
emitted photons has also been identified as a possibly 
efficient tool for discriminating between vacuum-effect- 
based models of SL and thermal light emission models. 

We argue, both here and elsewhere j||(J, that Casimir- 
like mechanisms for SL are viable, that they make both 
qualitative and quantitative predictions, and that they 
are now sufficiently well defined to be experimentally fal- 
sifiable. 
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